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We describe a measurement-based state preparation scheme for the efficient build up of cluster 
states in atom-cavity systems. As in a recent proposal for the generation of maximally entangled 
atom pairs [Metz et al, Phys. Rev. Lett. 97, 040503 (2006)], we use an electron shelving technique 
to avoid the necessity for the detection of single photons. Instead, the successful fusion of smaller 
into larger clusters is heralded by an easy-to-detect macroscopic fluorescence signal. High fidelities 
are achieved even in the vicinity of the bad cavity limit and are essentially independent of the 
concrete size of the system parameters. 
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I. INTRODUCTION 

In 2001 Raussendorf and Briegel pointed out that cer- 
tain highly entangled states present an innovative ap- 
proach to quantum computing |l| . The attractiveness 
of these so-called cluster states [3 arises from the fact 
that they can be grown off-line in a probabilistic fashion. 
Afterwards, a so-called one-way quantum computation 
can be carried out without having to create additional 
entanglement. Any quantum algorithm can then be per- 
formed using only single-qubit rotations and single-qubit 
measurements. Scalable fault-tolerant one-way compu- 
tation is possible, provided the noise in the implemen- 
tation is below a certain threshold [1, 1^. For exam- 
ple, Raussendorf et al. Q recently introduced a fault- 
tolerant three dimensional cluster state quantum com- 
puter based on methods of topological quantum error 
correction. Other authors identified highly efficient clus- 
ter state purification protocols [1, 0| . 

A very efficient way to create a cluster state of a very 
large number of atoms with very few steps is to employ 
cold controlled collisions within optical lattices with one 
atom on each site Q- Using this approach, Mandel et 
al. 0] already created cluster state entanglement and 
reported the observation of coherence of an atom de- 
localised over many sites. Unfortunately, single-qubit ro- 
tations cannot be easily realised, since laser fields applied 
to one atom generally affect also its neighbours. To fa- 
cilitate one-way quantum computing in optical lattices 
several schemes have been proposed for the realisation 
of single-qubit rotations without having to address the 
atoms individually [l^, [HI, . 

However, higher fidelities can be obtained using a 
measurement-based cluster state growth approach. An 
example is the linear optics proposal by Browne and 
Rudolph [l^. Using linear optics, a four-photon cluster 
state has already been generated in the laboratory p^ . 
Currently, the scalability of this approach is hampered 
by the lack of reliable photon storage. To overcome this 
and the above mentioned addressability problem in opti- 
cal lattices, quantum computing architectures have been 
proposed using hybrid systems based on atomic and pho- 
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FIG. 1: (Colour online) (a) Experimental setup for the cluster 
state growth of atomic qubits. Fusing two clusters requires 
placing an atom from each cluster into the resonator, where 
they are illuminated by a laser field, (b) Level scheme of 
a single atom in the cavity. Both atoms should experience 
similar coupling constants. 



tonic qubits [H, [H, [O, [H, [Til . To create entanglement 
between distant atoms, the atoms are operated as sources 
for the generation of single photons on demand followed 
by photon pair measurements in a carefully chosen basis. 
The main limitation of these two approaches lies in the 
difficulty of detecting single photons. 

Here we avoid the necessity of detecting single photons. 
Instead the successful fusion of two smaller clusters into 
a larger one is heralded by a macroscopic fluorescence 
signal. The experimental setup is shown in Fig. [ija). 
Suppose two atoms are placed into the antinodes of an 
optical cavity and a laser with Rabi frequency Q and de- 
tuning A is applied. It drives the 1-2 transition of each 
atom as shown in Fig. [Ijb) . The same transition should 
also interact with the resonator field with the atom-cavity 
coupling constant g and detuning A. We denote the de- 
cay rate of a single photon through the cavity mirrors by 
K [20| and F is the spontaneous decay rate of the excited 
atomic state. Here we are especially interested in the 
parameter regime where 



A > 17, g, K, F . 



(1) 



A detector constantly monitors the emission of photons 
through the cavity mirrors. 
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As in Refs. [2ll . |23 | , we assume that both atoms expe- 
rience similar interactions. Here the state |0) is decou- 
pled from any dynamics of the system. As we see below, 
there are therefore three distinct fluorescence levels in the 
leakage of photons through the cavity mirrors. These are 
similar to the discrete fluorescence levels of two dipole 
interacting three-level atoms, which exhibit macroscopic 
quantum jumps [23l. [2^. [25j . The lack of cavity photons 
indicates that both atoms are in |0). The emission of 
cavity photons at a lower intensity level indicates that 
one atom is in |0) and one atom is in |1) without reveal- 
ing which one, while cavity photons at a maximum rate 
indicate that the atoms are in |11). In fact, the observa- 
tion of the cavity fluorescence implements a probabilistic 
parity measurement with the projections 
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The successful projection of the atoms onto the subspace 
spanned by the states 101) and |10) can be used to create 
entanglement [1^, [13, [1^. For example. 



(Poi+Pio)[i(|0)-l-|l))®(|0) + |l))] 
= i(|01) + |10)). 



(3) 



However, it can also be used to generate entanglement 
between atoms without destroying any previous entan- 
glement of these atoms with other atoms. For example, 
the projection Pqi + Pio applied to atoms 2 and 3 ob- 
tained from two different Bell pairs. 
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results in the generation of a four-atom GHZ state. 
Browne and Rudolph moreover showed that the measure- 
ment ([2]) enables the fusion of two smaller cluster states 
into one larger one with a success rate of 50 % [l^ . It 
can therefore be used for the sequential build up of large 
cluster states. Detailed analyses on the scalability of re- 
lated probabilistic cluster state growth schemes can be 
found for example in Refs. [13, HS,, ,29,, ,3(1, ,3L] . 

Achieving high fidelities is possible, even when using 
moderate atom-cavity systems with relatively large spon- 
taneous decay rates. The reason is that the qubits are 
encoded in lon g-living atomic ground states. Moreover, 
as in Refs. [2l]. 12^. [271 f28j . we can allow for single atom- 
cooperativity parameters C, 



C's. For example, if 77 = 0.2 the cooperativity parameter 
C should be 5 or larger. 

The above described distinct fluorescence levels occur 
in the emission from the cavity mode for a very wide 
range of experimental parameters. The performance of 
the proposed state preparation scheme is therefore es- 
sentially independent of the concrete size of the system 
parameters. To illustrate this we show that high fi- 
delity parity measurements arc possible even when the 
two atoms experience coupling constants differ from each 
other by up to 30%. Once a cluster state has been 
built, performing a one-way quantum computation re- 
quires only single-qubit rotations and measurements as 
they are routinely used in ion trap experiments (ssl . Is^ . 
More specifically, read out measurements are performed 
via the creation of macroscopic fluorescence signals and 
too have a very high accuracy even when using finite ef- 
ficiency photon detectors [i^l . 

It is experimentally feasible to trap two atoms fairly ac- 
curately in different antinodes of the cavity field. The ef- 
ficiency of cavity cooling has recently been demonstrated 
by Nufimann et al. j41| . Domokos and Ritsch showed 
that it is possible to take advantage of cavity-mediated 
forces to keep the atoms predominantly at positions with 
maximum atom-cavity couplings [4^ . A disadvantage of 
the proposed state preparation scheme lies in the neces- 
sity to move atoms in an out of an optical cavity. The 
so-called shuttling of atoms [i^l is relatively time con- 
suming and limits the efficiency with which one can build 
large cluster states. However, its feasibility has already 
been demonstrated by sever al group s who combined for 
exainple atom trapping [s^, HI, or ion trapping 

[H, [43] technology with optical cavities. New perspec- 
tives arise from the development of atom-cavity systems 
mounted on atom chips [H, \^ ■ 

There are six sections in this paper. In Section [IT] we 
describe the setup shown in Fig. [T] and derive its effec- 
tive dynamics. In Section [Hi] we discuss the nature of 
the three distinct levels in the fluorescence through the 
cavity mirrors and describe how to exploit these for the 
implementation of the probabilistic parity measurement 
([2]) . In Section IIVI we analyse the performance of the 
proposed protocol. In Section |V| we review the cluster 
state build up with parity checks. Finally, we summarise 
our findings in Section IVII Some mathematical details 
are given in Appendix [^ 
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of the order of one and larger, as they are currently be- 
coming available in the laboratory [H, [H, [13, [H, [la, [13] ■ 
For C = 1 and when using a perfect single photon detec- 
tor, we show that it is possible to achieve fidelities above 
0.88. Lower photon detector efficiencies 77 require larger 



II. THEORETICAL MODEL 



In this section we use the quantum jump approach [49|, 
[sol . [51] ] to obtain an effective theoretical model for the 
description of the atom-cavity system in Fig. [TJ 
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A. The no-photon evolution 

Proceeding as in Ref. [49| . i.e. assuming rapidly re- 
peated environment-induced measurements and starting 
from the total Hamiltonian for the atom-cavity system 
and the surrounding free radiation fields, one can show 
that the (unnormalised) state of the system under the 
condition of no photon emission within (0, t) can be writ- 
ten as 



(6) 



Here |V'o) is the state of the system a.t t — and 
Ucond{t,0) is the no-photon evolution operator. The cor- 
responding conditional Hamiltonian equals, with respect 
to an appropriately chosen interaction picture. 
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The non-Hermitian terms in the last line of this equation 
damp away population in states that can cause an emis- 
sion. After renormalisation of the state vector \ip'^{t)), 
this results in a relative increase in population of states 
with a lower spontaneous decay rate. In this way, the 
quantum jump approach takes into account that the ob- 
servation of no photons reveals information about the 
system. It gradually reveals that the system is more likely 
to be in a state where it cannot emit. 

In the following, we decompose in Eq. 



as 



j,k=0 n=0 



(8) 



where Cjfc;„ is the amplitude of the state |jfc; n) with the 
first atom in \j), the second atom in \k) and n photons in 
the cavity mode. According to the Schrodinger equation, 
the evolution of these coefficients is given by 
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Writing out these equations and using Eq. ([T]) we find 
that the coefficients of states with population in |2) 
change on a much faster time scale than the coefficients of 
atomic ground states. We can therefore eliminate them 
adiabatically from the system's evolution. Doing so, and 
setting their derivative equal to zero, we obtain 
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up to second order in 1/A, given that most of the popu- 
lation remains in the atomic ground states. Substituting 
Eq. (fTO|l into the differential equation (l9|), we then find 
that 



COO;n — 

C01;n = 

ClO;ri = 

Cll;n = 
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up to second order in 1/A. Given the parameter regime different time scales. Since k is much larger than all other 
in Eq. ([T|), these differential equations contain two very frequencies that scale as 1/A or 1/A^, the coefficients of 
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states with n > 1 evolve much faster than the eoeffieients 
of states with n = 0. This aUows us to ehminate the cav- 
ity field adiabatically from the evolution of the system. 
Setting the derivative of the coefficients with n = I equal 
to zero, we find that 



C00;1 = 0, 



C01;l 



ClO;l 



Cll:l = 
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up to second order in 1/A. On average there is much less 
than one photon in the cavity mode. 

Finally we derive a set of differential equations for the 
ground state coefficients of the atom-cavity system. In- 
troducing the effective parameters 



^-^ff = iA 
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and substituting Eqs. pO)) and (fT2|) into Eq. ([9]), we ob- 
tain the effective differential equations 

C00;0 = 0, 

C01;0 = (iAoff — ^Foff — ^Koff) Coi;0 , 
ClO;0 = (iAoff — ^roff — 5«:cff) CiO;0 : 

cii;0 = (2iAcff - Tcff - 2Koff)cii;o • (14) 



Here the spontaneous decay rates are correct up to second 
order in 1/A, while level shifts small compared to Aos 
have been neglected. Eq. pi]) can be sumarised in the 
effective Hamiltonian 

iJeff = -^(Aeff + ireff + iKeff)[|01)(01| + |10)(10|] 

-ft(2Aeff + iEeff + 2iKeff) |11)(11| , (15) 

which acts only on the atomic ground states. The state 
1 00) is effectively decoupled from the dynamics of the 
system, while the states |01) and |10) evolve in the same 
way. Both cause an atomic emission with the effective 
decay rate Teg or the leakage of a photon through the 
cavity mirrors with n^s- The state |11) causes an atomic 
emission with the decay rate 2roff and the emission of a 
cavity photon with 4Koff. 



B. The effect of a photon emission 

The effect of a photon emission on the state of the 
atom-cavity system can be described with the help of 
supplementary jump or reset operators Rx- If \^) is the 
state prior to a photon emission of type x, then Rx\'4') 
is the (unnormalised) state immediately afterwards. For 
convenience we define the reset operators in the fol- 
lowing such that 



is the probability density for the corresponding emission 
to take place. 

According to the quantum jump approach [49j . the re- 
set operator for the emission of a photon via the 2-j 
transition of the atoms is given by 



R, 
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if the photons from the 2-0 and the 2-1 transition are dis- 
tinguishable. The discussion in the previous subsection 
shows that the only states with population in the excited 
atomic state are |02;0), |20;0), |12;0), |21;0), and |22;0). 
From Eq. (fTO)) we see that that coefficients of these states 
depend to first order in 1/A only on the coefficients coo;0, 
coi;0, cio;o and cn-o. Combining Eqs. (jlO|l and (fT7|) . we 
see that the reset operators for the photon emission from 
the atoms are effectively and up to an overall phase factor 
given by 



i?cff;o = v/T^[|00)(01| + |00)(10| + |01)(ll| 

+|io)(ii|], 

i?cff;i = [|01)(01| + |10)(10| + 2|11)(11|] (18) 

with 
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and Feffio + ^eS;i = Feff. Again we find that the states 
|01) and |10) have the atomic decay rate Fos, while |11) 
has the atomic decay rate 2Foff. 

During the leakage of a photon through the cavity mir- 
rors with decay rate k, one photon is removed from the 
resonator field. The corresponding reset operator is ther- 
fore given by 



i?c ~ ^fnh . 



(20) 
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In the previous section we have seen that there is on 
average at most one photon in the cavity. Only the states 
|01; 1), |01; 1), and |11; 1) contribute to a cavity photon 
emission. We know that their coefficients depend only on 
coi;Oi cio;o and cii;o. Combining Eqs. (fT^ and (|20p we 
find that the cavity jump operator Rq is to first order in 
1/A given by 

i?cff;c = \/;^;i[|01)(01| + |10)(10| + 2|ll)(ll|] . (21) 

Again, we see that the states |01) and |10) can cause a 
cavity photon emission with decay rate Koff, while |11) 
has the collectively enhanced cavity decay rate 4Keff. 



C. The master equation 

It should also be noted that the quantum jump ap- 
proach above is consistent with the master equation, 
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FIG. 2: Illustration of the effective evolution of the atoms. 
When spontaneous emission from the atoms is negligible, the 
state of the atoms converges within a relatively short time into 
one of the three subspaces D, L, and H . Each subspace is a 
characterised by a discrete fluorescence level in the leakage of 
photons through the cavity mirrors. 



which is often alternatively used for the description of 
an open quantum system. It reads 



P 



7^cfr(p), 



(22) 



and can be obtained by averaging over all the possible 
traj ectories that a single atom-cavity system can undergo 
[49|. The master equation is particularly well suited to 
the prediction of ensemble averages. For example, in- 
stead of displaying the existence of discrete levels in the 
fluorescence of a single system, the master equation can 
be used to predict the intensity of the emitted light av- 
eraged over all possible trajectories. 



III. BASIC PRINCIPLE 

Let us now discuss the dynamics of the atom-cavity 
system shown in Fig.[T]in more detail. Suppose the atoms 
are initially in 

\i'o) = coo|00)+coi|01)+cio|10)+cn|ll). (23) 

As we see below, there are three distinct fluorescence 
levels in the leakage of photons through the cavity mir- 
rors. Their origin is the existence of three decoupled 
subspaces in the effective evolution of the atomic ground 
states, when spontaneous emission from the atoms re- 
mains negligible. We denote them in the following by 
D, L and H (c.f. Fig. [J). The emission of photons at a 
certain rate or their complete absence gradually reveals 
information about the atoms. This gain of information 
gradually increases the population in one subspace with 
respect to the others until the population in any other 
subspace is irreversibly lost. The result is the projection 
of the atomic state IiJjq) into one subspace. 



FIG. 3: (Colour online) A possible trajectory of the atom- 
cavity system given the initial state Q obtained from a quan- 
tum jump simulation with To = Ti = 0.1 k, g — Q = it, and 
A = 50 K. Here no photons are emitted and all population 
eventually accumulates in the D-subspace. 



no cavity photon is emitted for a time of the order of 
a few 1/Koff, we therefore learn that the atoms are in 
a state where they cannot transfer population into the 
cavity mode. Consequently, the relative population in 
1 00) increases, while the population in the L and in the 
_ff-subspace decrease. Eventually this results in the pro- 
jection of the state ((^5]) into |00) as shown in Fig. [31 

The occurrence of two distinct fluorescence periods 
when the atoms do lead to the emission of cavity photons 
may seem less obvious. Figs. S] and [S] show possible tra- 
jectories of the system in these cases. Any population in 
1 00) vanishes with the first photon emission. Fig. 2] shows 
how relatively frequent cavity emissions result in an ever 
increasing population in |11), since this reveals that the 
atoms are most likely in the ff-subspace. Eventually all 
atomic population accumulates in |11), and can no longer 
return into another subspace via cavity photon emission. 
The system therefore continues to emit photons at its 
maximum rate given by 4Kcff. 

However, it is also possible that the time between two 
subsequent photons is comparatively long. Such less fre- 
quent events result in a relative increase in population of 
the states |01) and |10) with respect to the population in 
1 11). This is illustrated in Fig. [S] The reason for this is 
that seeing no photon for a time that is long compared 
to 1 /4Koff and after the first photon emission reveals that 
the system is more likely to be in the subspace with the 
lower emission rate. Eventually, the population in |11) 
vanishes completely and the system emits cavity pho- 
tons at the rate k^s- Numerical simulations confirm that 
the probability of the odd-parity projection, Pqi + Pio 
(c.f. Eq. (21)), indeed equals the initial population in the 
L-subspace, as predicted for an ideal measurement. 



B. Implementation of a parity measurement 



A. Three distinct fluorescence levels 

Here the laser interaction is chosen such that the state 
1 00) is not involved in the evolution of the system. If 



In order to implement the parity measurement ^ in 
the absence of spontaneous emission from the atoms, it 
is sufficient to observe the cavity fluorescence over a time 
T long enough to clearly distinguish the three fluores- 
cence levels mentioned above. The concrete protocol is 
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FIG. 4: (Colour online) A possible trajectory of the atom- 
cavity system for the same initial state and the same param- 
eters as in Fig. [S] Here photons are emitted at a relatively 
high rate. This results eventually in a projection of the atomic 
state into the i^-subspace. 
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FIG. 5: (Colour online) A possible trajectory of the atom- 
cavity system for the same initial state and the same param- 
eters as in Figs. [3] and |4] Here the time between two photon 
emissions is initially relatively long. This reveals that the 
atoms are most likely in the L-subspace and all population in 
D and H gradually vanishes. 



as follows: 

1. Place respective atoms inside the resonator. 

2. Turn on driving laser for a time T longer than a few 
1/Koff and count the number of detected photons. 

The parity measurement is successful and results in the 
projection Pqi-I-Pio, when the number of emitted photons 
is close to rjKosT, which is the average number of photon 
detections when the system is in the L-subspace. The 
average number of photon detections when the system 
is in the i/-subspace is given by ArjKcsT. No photons 
indicate a projection onto the D-subspace. Fig. [5] shows 
the fidelity and event probability for different detector 
click events for operation ([3]) and for a concrete choice of 
experimental parameters (C = 10). Even in the presence 
of a non- negligible spontaneous decay rate of the atoms, 
it is possible to achieve fidelities well above 0.9. 
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FIG. 6: (Colour online) Fidelity of the prepared state in case 
of n detector clicks in (0, T) averaged over many trajectories 
and the corresponding probability for this to happen. Here we 
compare the results from an analytical calculation [s^j with 
the results obtained from a quantum trajectory simulation 
with T = 5/Kcff, Fq = Pi = 0.05k, g = Q = k, A = 50k, and 
77 = 1. A photon number close the KcffT indicates that the 
atoms are most likely in the L subspace, while no detector 
clicks or a relatively large number of clicks shows that the 
atoms are in |00) or respectively. 



C. Optimised protocol 

Deviations in the fidelity from unity in Fig. [S] are 
largely due to the possibility of photon emission from the 
atoms with Feff. For example, when starting in |11), the 
emission of a photon with reset operator i?eff;0 in Eq. (flSl) 
projects the atoms onto the L-subspace. Similarly, such 
an emission can transfer population from the L-subspace 
into the D-subspace. In both cases, it might be assumed 
that the atoms are in L, consequently resulting in a de- 
crease of the fidelity of the prepared state. When using 
the parity measurement for the generation of multi-qubit 
entanglement, as described for example in Eq. (|3|), the 
result is the loss of the entanglement with atoms outside 
the cavity. 

Reducing the occurrence of atomic emissions requires 
shortening the interaction time T. However, in the above 
protocol, this would make it difficult to distinguish the 
fluorescence when in H from that in L. As a solution we 
propose to use the double heralding technique of Barrett 
and Kok We now consider the following variant of 
the above protocol (c.f. Fig. [7]): 

1. Place respective atoms inside the resonator. 



2. Turn on driving laser for a maximum time 
until the first detection of a photon. 

3. Swap the states |0) and |1) of each atom. 



or 
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FIG. 7: (Colour online) Schematic view of the main stages 
of the optimised protocol. After the atoms have been placed 
into the resonator in step 1, the laser field is switched on in 
step 2 until a cavity photon is detected at time T\. In step 
3 a TT-puIse transfers |0) into |1) and vice- versa. In step 4 
the laser is turned on again until a second click at T2. This 
second click signals an odd parity state of the atoms. 



4. Repeat step 2 and 3. 

This protocol allows us to measure the parity of atoms 
in a relatively short time. In the ideal case, the only 
events that produce two detector clicks are due to the 
system being in L. The reason for this is that atoms in 
1 00) cannot produce a click in step 2. Atoms initially in 
1 11) are transferred into |00) in step 3 and are therefore 
unable to emit a photon in step 4. 

Finally we comment on the optimal size of Tmax- The 
duration of the applied laser pulse should be comparable 
to the mean time between two photon detections when 
the system is in the L-subspace. This is sufficient to as- 
sure that there is a high probability for the detection of a 
photon in step 2 and 4, if the atoms are in L. For exam- 
ple, if Tniax = 3/?7Kcff7 this probability is already above 
90 %. Longer laser pulses do not increase this probability 
significantly and can therefore be avoided. They might 
only lead to a slight decrease of the fidelity of the odd- 
parity measurement due to an increase in the probability 
for the spontaneous emission of a photon from one of the 
atoms before the detection of the first cavity photons. 
In the remainder of this paper, we nevertheless assume 
that Tmax ~ oo. This simplifies the following calculations 
which nevertheless yield good approximations for the ac- 
tual fidelity and the success rate of the proposed scheme 
for finite T,„„^. 



IV. PERFORMANCE ANALYSIS OF THE 
OPTIMISED PROTOCOL 

In the following we analyse the optimised protocol in 
detail and show that its performance is comparable to the 
scheme presented in Ref. [2l| . As an example, we con- 
sider the operation described in Eq. Q and calculate the 
average fidelity and success rate for the preparation of a 
four-atom GHZ state. In Section IIV Al we analyse the 
ideal scenario of perfect photon detection. Finite pho- 
ton detector efficiencies are taken into account in Section 



IIV Bl In Section llV Cl we emphasize that the performance 
of the proposed state preparation scheme is essentially in- 
dependent of the concrete size of the system parameters 
and hence very robust against parameter fluctuations. 



A. Average fidelity for unit efficiency photon 
detectors 

One factor that decreases the fidelity of the prepared 
state is population in excited atomic states after the laser 
field has been turned off. This population might result 
in an atomic emission, which transfers the system into a 
state with a reduced overlap with the target state. From 
Eq. (fTOj) we see that it equals 



atom excited 



4A2 



(24) 



up to second order in 1/A, in the case of the odd par- 
ity projection Pqi -I- Pio. This population can be made 
arbitrarily small even in the presence of relatively large 
spontaneous decay rates by simply increasing the detun- 
ing A. In the parameter regime ([T]), corrections due to 
the population in Eq. are hence negligible. Possi- 
ble remaining cavity excitations do not affect the fidelity 
of operation ^ , since their emission does not affect the 
state of the atoms once the projection onto one of the 
subspaces D, L ot H has occurred. 

Corrections to the fidelity of the state prepared 
through an odd-parity projection are generally domi- 
nated by effects due to spontaneous emission from the 
atoms. To calculate these corrections for unit efficiency 
photon detectors, we consider the following two events: 

• Event A: The atoms emit the first cavity photon 
at Ti in step 2 and at T2 in step 4. Spontaneous 
emissions from excited atomic states may occur but 
the atoms are finally in the desired state. 

• Event B: The atoms emit the first cavity photon 
at Ti in step 2 and at T2 in step 4. Spontaneous 
emissions from excited atomic states occur and the 
final state of the system is not the desired one. 

If wx{Ti,T2) denotes the probability density for having 
event X, the fidelity of a successful odd parity check 
equals 



FiTi,T2) 



WAiTl,T2) 



Wa{Ti,T2) +Wb{Ti,T2) 



(25) 



We now analyse this fidelity in terms of the probabilities 
of Markovian processes. 

To do so, we denote the probability density for the 
spontaneous emission of a photon from the atoms, which 
transfers the atoms from the X-subspace into the Y- 
subspace, by ^xy- Moreover, kx is the probability den- 
sity for a cavity photon emission, when the atoms are 
in X. From Eq. (PT|) we see that such an emission does 
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not change the state of the atoms within the respective 
subspace. Given an initial state, where the coefficients 
of the |00), |01), |fO) and components of the state 
of the atoms inside the cavity are initially all the same 
(as in Eq. (|4])), the probabihty density wa{Ti,T2) can be 
written as 



oo 

X ^ e-"^^^ pL^HT2) F{n + m) 



m=0 



n=0 



(26) 



Here the factor ^ is the initial population in the L- 
subspace and e~'^^'^ is the probability for no cavity de- 
cay in (0, T) given that the atoms are in L. Moreover, 

Pn"\T) is the probability of n atomic emissions from 
subspace L back into L and no emissions from L into 
D in (0,T). An expression for this probability can be 
found in Eq. (|A2p . The fidelity of the final state after 
n such spontaneous emissions is denoted by F{n). Here 
F{0) = 1 and F{n > 0) = ^ due to the nature of spon- 
taneous emission from the atoms. 

Similarly we now calculate the probability density for 
event B to occur. Doing so we find that 



wb{Ti,T2) = \ dt E ^HL e-«-*Pif ) it) 

oo 

X E 4 g-K,(Ti-t-hT.) piL)^rj.^ -t + T2). 



n=0 



+-2 E e---(^^+^^) )(Ti + r2)(l - Fin)) . 

(27) 



ri=0 



The factor j in the first line is the initial population in the 
-ff-subspace. Here t is the time at which a transition from 
H to L occurs. To calculate wb(T'i, T2), we integrate 
over all possible values for t. Moreover, Q-hHL+KH)t 
the probability for no atomic emission out of H and no 
cavity decay in (0, t) given that the atoms are in H. In 

analogy to the above notation, Pn ' (t) is the probability 
of n atomic emissions from subspace H back into H and 
no emissions from H into L in (0,i). An expression for 
this probability can be found in Eq. (|A3[) . The last line 
in Eq. ((27|) is analoguos to Eq. ((26)) and contains the 
probability density that the final state of the atoms does 
not overlap with the desired state although the atoms 
were initially in L in the absence of a transition from the 
H into the L-subspace. 

Using Eqs. (|Aip and (|A3p . the probability densities 
wa{Ti,T2) and wb{Ti, T2) can be calculated analytically. 



Doing so. we obtain 



wa{Ti,T2) = jk\ 



^-(hl+1ld+1ll){Ti+T.2) 



and 

Wb{Ti,T2) 



IHL 

Kh - Kl + -fHL - ILD 



+T2) _ q~(k.l+1Ld)T2 ^-{KH+'rHL)Tj 



(28) 



■■) 

(29) 



It is possible to optimise the corresponding fidelity 
F{Ti,T2) in Eq. (^5]) by postselccting events, where Ti 
and T2 are both short. In such cases, the probability for 
a decrease of the fidelity due to an atomic emission is low 
and wb{Ti, T2) remains negligible. However, this optimi- 
sation comes at the cost of a decrease of the success rate 
for an odd parity check. 

To maximise the efficiency of the proposed cluster state 
preparation scheme, let us accept all events independent 
of the size Ti and T2 and assume that the interaction time 
Tmax in Fig. [7] is very large. The average fidelity of the 
final state in the case of a projection onto the L-subspace 
is then given by 



Pa- 



Pa 



Pa + Pb 



(30) 



in analogy to Eq. (pS)) . if Px is the probability for an 
event X = A, B to take place. Moreover 



Psuc = Pa + Pe 



(31) 



is the probability for the observation of such a projection 
and equals the denominator in Eq. (j30p . Since 



Px — lim 
we find, using Eqs 
Pa 
Pb 



dPi / dT2Wx{Ti,T2), (32) 







4^L 



and 
1 



1 



IHL 



iLhilLL + 2kl + 27LI3) 



(kl + ILD + ILl)"^ 



KH + IHL _ 

with the transition rates (c.f. Section [H]) 

IHL = '^ILL = 2reft;l , 

IHH = 27l£) ~ 2roff;0 , 



(33) 



(34) 
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FIG. 8: (Colour online) Logarithmic plot of the fidelity of 
the output state against the detector efficiency 77 for several 
values of the cooperativity parameter C. The lines show the 
analytical solution given in Eq. (|38p for the four-atom GHZ 
state preparation described in Eq. Q, while the symbols have 
been obtained from a quantum trajectory simulation of the 
simpler entangling operation (|3} assuming A = 50 k, and g = 
Q = K. The triangle, circle , square, and diamond, correspond 
to the C = 5, 10, 20, and 40, respectively. 



In the special case, where Tq 
pUfy to 



Fi, these equations sim- 



27LL = 27l£1 = "fHL = IHH , 

KH = 4kl = 16C -^HH , (35) 

and the average fidelity and the success rate depend only 
on the single atom cooperativity parameter C, 



and 



5/32 + 4C + 28C2 + 64C3 
3/8 + 7C + 38C2 + 64C3 



6C^ + 64C3 
1/8 + 4C + 40C2 + 128C3 



(36) 



(37) 



In the limit of large C, the fidelity i^av of the entan- 
gling operation ^ approaches unity, while the success 
rate Pguc converges to ^, as expected for an ideal pro- 
jection. In the presence of non-negligible atomic emis- 
sion, i.e. for smaller C's, the fidelity and success rate are 
slightly smaller. Even when C = 1, wc have Fav = 0.88 
and Psuc = 0.41 in the presence of ideal photon detectors 
with 77 = 1. Average fidelities exceeding 0.99 become 
possible when C approaches 20. 



Average fidelity for finite efficiency plioton 
detectors 



To calculate the average fidelity and success rate for 
an odd-parity projection for finite photon detector effi- 
ciencies ?7 < 1, we consider the events: 
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FIG. 9: (Colour online) Probability of a successful odd parity 
check as a function of 77, for several values of C and for the 
same parameters as in Fig. [S] The circle , square, cross, and 
triangle, correspond to the C = 5, 10, 20, and 40, respectively. 
The lines show the analytical solution given in Eq. (|39|) for 
the four-atom GHZ state preparation described in Eq. Q. 



• Event A': The first detection of a cavity photon 
at Ti in step 2 and at T2 in step 4. Spontaneous 
emissions from excited atomic states may occur but 
the atoms arc finally in the desired state. 

• Event B': The first detection of a cavity photon 
at Ti in step 2 and at T2 in step 4. Spontaneous 
emissions from excited atomic states occur and the 
final state of the system is not the desired one. 

We then notice that rj has no effect on the probabil- 
ity densities for atomic emissions. It only extends the 
mean time until the detection of a first cavity pho- 
tons. To obtain the probability densities wa'{Ti,T2) 
and wb'{Ti,T2), we therefore only need to replace kx 
in Eqs. ^ and ^ by rjKx and e-*"^* by e"''''^*. Pro- 
ceeding as above, we then find that the average fidelity 
and success rate are now given by 



F — 

a,v 



and 



5/32 + A-qC + 28(77C)2 + 6A{i]C)^ 
3/8 + 7r]C + 38(?7C)2 + 64{t]C)^ 



6{r]C)^ + 64(??C)3 
1/8 + 4j]C + 40(?7C)2 + 128(77(7)3 



(38) 



(39) 



These arc the same expressions as in Eqs. (|36l) and ([57)) 
but with C now replaced by 77c. 

Figs. [5] and [H] compare the average fidelity in Eq. (pS)) 
and the success rate for an odd-parity projection in 
Eq. p9p with the results obtained from a numerical so- 
lution of the evolution of the system using Eqs. ([7]), 
(fT7)) . and (PO]) for the entangling operation described in 
Eq. ([3]). Operation ([3]) is slightly different from opera- 
tion which wc considered in the above calculations, 
since atomic emissions lead less frequently to an error in 
the simpler pair entangling scheme, as there is no initial 
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entanglement that needs to be preserved. Nevertheless, 
there is relatively good agreement between both curves. 
Figs. [8] and [9] indeed confirm that reducing rj has the same 
effect as replacing C by 77C. 



C. Parameter dependence 

Due to its postselective nature, the performance of the 
proposed state preparation scheme is essentially indepen- 
dent of the concrete system parameters. Fig.[8]shows that 
fidelities well above 99 % are inevitably when rjC > 20. 
Whenever this condition is fulfilled, there are three dis- 
tinct fluorescence levels in the emission of cavity photons. 
The scheme is constructed such that the medium level al- 
ways indicates that one atom is in |0) and one atom is 
in |1) without revealing which one. Turning off the ap- 
plied laser field upon the detection of a photon in step 2 
and 4 is hence sufficient to realise the parity operation ([2]) 
with very high accuracy. The proposed state preparation 
scheme is therefore robust against a parameter fluctua- 
tions, like moderate fluctuations of cavity coupling con- 
stants and laser Rabi frequencies. 

In the following, we show that high fidelities are 
achieved even when the atoms experience quite differ- 
ent coupling constants. As an example, we consider the 
case where the cavity-coupling constant of atom 1 is gi 
and the cavity-coupling constant of atoms 2 is given by 
g2- For simplicity we neglect spontaneous emission from 
the atoms (F = 0) in the following. Proceeding as in 
Section ini we find that the conditional Hamiltonian ([T5|) 
is now given by 



eft = -;i(Aeff +^Kcff;2) |01)(01| 
-;i(Aeff + iKcff;l) |10)(10| 

-;i[2Aeff + i(V7^;iT+ |11)(11| (40) 

with 

- (41) 

At the same time, the reset operator in Eq. (|2T|) becomes 

Re«-c = V7J^(|10)(10| + |11)(11|) 

+^/^(|01)(01| + |ll)(ll|). (42) 

These two equations can be used to simulate for example 
all the possible trajectories for the two-qubit entangling 
operation ([3]). 

Let us first consider the case, where the atoms perma- 
nently see different coupling constants gi. Fig. 1101 shows 
the average fidelity F^v of the prepared state in the case 
of an odd parity projection as a function of e, where 



KcS-l = (1 + e) Koff , KcS-2 = (1 - e) Kcff : 



(43) 



Fidelities above 0.9 are achievable, even when Keff;i and 
Keff;2 differ by their average value Koff! The reason for 
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FIG. 10; (Colour online) Average fidelity, Fav, from Eq. (|49)) 
for the entangling operation described in Eq. Q as a function 
of e. Fidelities above 0.9 are achievable, even when fi:cff;i and 
f^c«;2 differ from each other by Kcff. 



this is the postselective nature of the proposed optimised 
protocol. 

To illustrate this we now examine the dynamics of the 
system in the absence of spontaneous emissions from the 
atoms in more detail. The state of the atoms in case of 
a detector click at Ti and one at T2 is given by 



ReS;C C/cond(T2, 0) i?eff;C ?7cond (Ti , 0) \z^o) 



(44) 



with J/tt being the operation that exchanges the states 
|0) and |1). Using Eqs. (gHl) and gH), we find that the 
system is finally in the state 



l^(ri,T2)) 



1 



-(Koff;2Tl+Koff;lT2)/2 



101) 



Calculating the overlap of this state with the Bell state 
in Eq. ([3]), we obtain the fidelity 



F{Ti,T2) 



-(Krff;l-|-Kotf;2)(Tl-|-T2)/2 



2 Q— (Kctt;lTi+Kott;2T'2) _|_ g— ( Kotf ;2 T^i +reotf ; 1 12 ) 



(46) 



To calculate the average fidelity for the entangling op- 
eration ([U, we notice that the probability density for a 
click at Ti and T2 is in this case given by 



7(Tl,r2) = iKcff;lKcff;2 ( e-^''^^'^^^^ 



- ( Kcff ; 1 + Koff ; 2 T2 ) 



(47) 



The average fidelity i^av is again obtained by integrating 
over all possible click times Ti and T2, 

p _ S^dT,S^dT2w{T,,T2)F{T,,T2) 

J-dT,J-dT2wiT„T2) ■ ^''^ 
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FIG. 11: Fusion of two linear cluster states. A successful 
parity check creates a double-encoded qubit which links the 
two chains. If the parity check fails we project either on |00) 
or on 1 11). In both cases one qubit from each chain decouples. 



Inserting Eqs. ((43)) . (|46)) and (|47|) into this equation, we 
find that 



F - 1 - 

-£ av — ^ o ^ ■ 



(49) 



This means that very large fidelities are possible even for 
non-negligible e (c.f. Fig. [TO]). 

The very high fidelities in Fig. [TO] are due to the con- 
crete form of the state of the atoms after two photon 
emissions in Eq. ([^5)1 . As one can see, the states |01) 
and 1 10) have approximately the same coefficients, when 
Keff;2Ti -I- Keff;iT'2 and Kcs-iTi + Kes-2T2 are of about the 
same size. This applies for a very wide range of click 
times Ti and T2 and decay rates k^s-,! and KeS;2- Proceed- 
ing analogously, one can show that the proposed realisa- 
tion of the odd-parity check ^ is robust against temporal 
fluctuations of the atom-cavity coupling constant g and 
the laser Rabi frequency Q. 



V. CLUSTER STATE GROWTH 

In this section, we describe how to use the probabilistic 
parity check in Eq. ([2]) for the build up of two-dimensional 
cluster states 0- These highly entangled states consti- 
tute the main resource for one-way quantum computing 
[U, Once a cluster state has been built, local oper- 
ations and single-qubit read out measurements are suf- 
ficient to realise any possible quantum algorithm. That 
the projection ^ can be used to build cluster states has 
already been noted by Browne and Rudolph [l3|. Below 
we proceed in a similar fashion. 



A. Fusion of one-dimensional clusters 

Let us first have a closer look at the build up of a lin- 
ear cluster state, which consists of a chain of atoms with 
next neighbor entanglement. Larger clusters can be ob- 
tained through the fusion of two cluster states [l3| . In 
the scheme proposed here this requires placing one end 
atom from each chain into the optical cavity and perform- 
ing the parity measurement Entanglement between 
these two atoms is obtained and the original correlations 
with the atoms outside the cavity are preserved in case of 



a projection with Pqi + Piq. In case of a projection with 
Poo or Pii, the atoms in the cavity decouple from the rest 
and a new attempt has to be made to incorporate them 
into a larger cluster. 

More concretely, we now consider the case of fusing 
two cluster chains with n — m and m qubits, respectively, 
into one linear cluster of size n — 1. The initial state of 
the atoms can then be written as 



1 



2"/2 



i—7n~\-l 

(|0),+ari)|l),) 



Here tri*^ is the Pauli matrix 



(50) 



(51) 



with the exception of cri™' and ai^'' . These are given by 



a('")=af ^|0)(0| + |1)(1|. 



(52) 



By detecting odd parity for qubit m + 1 and qubit m we 
project these two qubits with Pqi + Pio. The resulting 
state is 



1 



2(n-l)/2 



) (|o), + 4-^^11)0 

i—ni-\-2 

)((7(™-i'|0)„+i|l)„ + |l)™+i|0). 



(53) 



The two originally independent chains are now linked 
via a double-encoded qubit, as illustrated in Fig. 1111 In 
order to remove the superfluous qubit we perform the 
Hadamard gate 



H 



1 

V2 



[|0)(0| + |0)(1| + |1)(0|-|1)(1|] (54) 



on atom m and measure its state. If qubit m is found 
in |0), we need to apply a operation to qubit m — 1 
to conclude the fusion of the two states. If we find qubit 
m in |1), we instead perform the cr^ operation on qubit 
TO — 1 and qubit to -|- 1. In both cases we end up in the 
n — 1 qubit cluster state 



IV') 



2(n-l)/2 



.(^-1)1 



(55) 



1=1 



Here we changed the indices of the atoms in order to 
close the gap caused by the loss of atom to. If instead 
the parity-check fails, and the projection Pqo or Pn is 
performed, then the qubits m + I and to are both pro- 
jected either on state |0) or |1). They are then decoupled 
from their respective cluster chains, which are now of 
length n — TO — 1 and to — 1. In order to increase the 
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FIG. 12: Fusion of two linear clusters into one two- 
dimensional cluster state. In case of a failure of the parity 
check, both chains split into shorter ones. 



efficiency of the growth of multi-qubit cluster states, one 
can abandon the measurement of atom m, as pointed out 
by Nielsen in Ref. [1^. Its presence can be used later to 
increase the success rate for a later fusion of a cluster to 
this part of the chain. 



B. Fusion of two-dimensional clusters 

Similarly, large two-dimensional cluster states are ob- 
tained via the fusion of two smaller structures into one. 
As a simple example, we now discuss the fusion of two 
one-dimensional cluster chains of length to and n into a 
simple two dimensional structure, as illustrated in Fig. 
W2\ The initial state of the two chains is given by 



1^) 



1 



2JV/2 



(g)(|0).+ari)|l),)(g)(|0),+ai^"-^)|l),) 
)),+a('=-i)|l),)...(|0),+ari)|l),)... 

(56) 



where TV = n + m is the total number of qubits. For 
example, if the parity projection Pqi +Pio is successfully 
applied to qubit fc and qubit Z, then these two qubits be- 
come a double encoded qubit which links the two chains. 
As a final step, one of the two atoms in the link, i.e. atom 
k or atom I, may be removed. In analogy to Section IVAl 
this can be done using again a Hadamard operation and 
a qubit read out measurement followed by operations. 

In the case of a Pqo or a Pn projection, qubit k 
and qubit I decouple from the rest of the cluster states, 
since their state is now known. The original chains 
split and the fusion of the two chains failed. The sit- 
uation is now worse than before. Instead of one large 
cluster we obtained four smaller ones and two single 
qubits. Nevertheless, it is possible to grow cluster states 
of any size with the help of probabilistic measurements. 
This applies even when the probability for the success- 



ful fusion of two clusters is below ^ [Ij]. More details 
about the scaling behaviour of similar probabilistic clus- 



ter state growth schemes can be found for example in 
Refs. [i3,H[2i,[33,[3l. 



VI. CONCLUSIONS 

In this paper we describe a scheme for the sequen- 
tial build up of atomic cluster states with the help of 
the probabilistic parity measurement ([2]). This measure- 
ment can be implemented via the detection of a macro- 
scopic fluorescence signal. It requires placing two atoms 
simultaneously into an optical cavity, where both ex- 
perience comparable cavity-coupling constants and con- 
stant laser driving with comparable Rabi frequencies 
(c.f. Fig. [T]). Fluorescence at a maximum level indicates 
that the atoms arc in |11), while fluorescence at a rel- 
atively low level indicates that the atoms are in |01) or 
1 10), without revealing which atom is in which state. In 
the case of no cavity photon emissions, the atoms project 
into 1 00). In Section [1111 we showed that the origin of 
these three distinct fluorescence levels is the existence of 
approximately decoupled subspaces in the effective evo- 
lution of the atomic ground states. 

One way to perform the parity measurement ([2) is to 
turn on the laser field for a fixed time T and to count 
the number of cavity photon detections in (0,T). How- 
ever, higher fidelities are achieved, when minimising the 
time T for which the laser field is turned on. This min- 
imises the effect of spontaneous emission from the atoms, 
which might disrupt the coherence between |01) and |10) 
or transfer population from |11) into a state with one 
atom in |0). Wc therefore propose an optimised protocol 
in Section lTlI C[ which makes use of the double heralding 
technique of Barrett and Kok [13]. It requires to turn 
off the laser field upon the detection of the first photon 
and to swap of the states |0) and |1) in both atoms. Af- 
terwards, another laser pulse is applied for a maximum 
time Tinax or until the detection of a second photon. In 
this way it is possible to measure how many atoms are 
in |0) in a much shorter time than in the first mentioned 
protocol. 

A detailed performance analysis of the optimised pro- 
tocol can be found in Section IIVI The main motivation 
for the proposed state preparation scheme is to allow for 
relatively large spontaneous decay rates and finite photon 
detector efficiencies 77. Indeed, it is possible to achieve fi- 
delities well above 0.99, when r\C > 20, while 77C > 1 is 
sufficient for fidelities above 0.88 (c.f. Fig. [5]). The suc- 
cess rate for an odd parity check is close to ^ for most 
detector efficiencies rj and values of the single atom co- 
operativity parameter C (c.f. Fig. [9]). This means, the 
performance of the proposed state preparation scheme 
is essentially independent of the concrete size of the ex- 
perimental parameters. Consequently the scheme is very 
robust against parameter fluctuations. To illustrate this, 
we show that the fidelity reduces by only 0.1 even when 
the effective atom-cavity coupling strengths both differ 
by approximately 30 % from their mean value. Fidelities 
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in excess of 0.99 of the values calculated for equal cou- 
pling constants require that the atom-cavity couplings 
differ by less than 10% (c.f. Fig.[10l). 

In Section |Vl we show how the parity measure- 
ment ([2|) can be used to grow two-dimensional cluster 
states. It has already been shown in the literature 
(c.f. e.g. Refs. [13, M, M,M,^) that the build up of 
large cluster states is possible even when the probability 
for the successful fusion of two clusters is below i . Here 
we propose a scheme, in which the success rate for an odd 
parity check is close to ^ even in the presence of finite 



efficiency photon detectors. Our cluster state growth 
scheme with macroscopic heralding is therefore expected 
to be much more practical than recent schemes based on 
the detection of single photons [H, [H, [O, [H, [H ■ 
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is in this subspace at t = and remains there throughout. 
It is given by 

JQ 

X / " dtn-1 ILL C-(^"+T-°)(*-*-i) ... 

Jo 

X ( ' dti -fLL e-(T^^+''^")*i , (Al) 
Jo 

if the ti denote the corresponding jump times. The evahi- 
ation of the above integrals is straightforward and yields 

Similarly, the probability for n atomic emissions from 
the H to the H subspace, given that the system is in this 
subspace at i = and remains there throughout, is given 
by 



APPENDIX A: CALCULATION OF Pi^'(t) AND Pn"\i) = / dt„... / dti-i^H^'^^""^'^" 

Pi."\t) ^° ^° 



a~(lHH+lHL) 



We now calculate the probability for n atomic emis- 
sions from the L to the L subspace, given that the system 



n 



